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ABSTRACT 

We are undertaking an imaging study of local star-forming galaxies in the Lyman-alpha (Lya) 
emission line using the Solar Blind Channel (SBC) of the Advanced Camera for Surveys onboard 
the Hubble Space Telescope. Observations have been obtained in Lya and H-alpha ( Ha) and six 



line- f ree continuum filters between ~ 1500A and the /—band. In a previous article ( Haves et al 
2005) we demonstrated that the production of Lya line-only images (i.e. continuum subtrac- 
tion) in the S^C-only data-set is non-trivial and that supporting data is a requirement. We 
here develop various methods of continuum subtraction and assess their relative performance 
using a variety of spectral energy distributions (SED) as input. We conclude that simple as- 
sumptions about the behavior of the ultraviolet continuum consistently lead to results that are 
wildly erroneous, and determine that a spectral fitting approach is essential. Moreover, fitting 
of a single component stellar or stellar+nebular spectrum is not always sufficient for realistic 
template SEDs and, in order to successfully recover the input observables, care must be taken 
to control the contribution of nebular gas and any underlying stellar population. Independent 
measurements of the metallicity must first be obtained, while details of the initial mass function 
play only a small role. We identify the need to bin together pixels in our data to obtain signal- 
to-noise (S/N) ratios of around 10 in each band before processing. At S/N = 10 we are able to 
recover Lya fluxes accurate to within around 30% for Lya lines with intrinsic equivalent width 
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(Why a ) of 10A. This accuracy improves to < 10% for VFL ya =100A. We describe the method of 
image processing applied to the observations presented in Ostlin et al. (2009) and the associated 
data-release. We also present simulations for an observing strategy for an alternative low-redshift 
Lya imaging campaign using ACS/SBC using adjacent combinations of long-pass filters to target 
slightly higher redshift. 

Subject headings: methods: data analysis — techniques: image processing — techniques: photometric 
— galaxies: starburst 



1. Introduction 

The Lyman-alpha emission line (Lya) is a pow- 
erful and frequently exploited observational signa- 
ture through which the galaxy population can be 
probed at high-redshift (z). In principle Lya can 
be used to probe the ionizatio n fraction during the 
final stages of re-ionizatio n (jMalhotra Sz Rhoadd 
12004 iDiikstra et all [2007h . cosmic star-formation 



rates (|Hu et al.ll998tlKudritzki et alj2000HAiiki et al 



2003) , large scal e stru cture (jVenemans et al 



20021: lOuchi et all 120051 ), and to identify po- 



tential host of populatio n ill star formatio n 
(jMalhotra fc Rhoadd 120021 : iNagao et all \200ii) . 
Ultimately, exactly how the Lya-emitting (or non- 
emitting) high-z galaxy population relates to the 
ultraviolet-selected Lyman break galaxy (LBG) 
population is uncertain, and studies of the high-z 
population of Lya emitters (LAEs) are interesting 
in their own right. 

In all surveys for which the ultimate science 
goal is more fundamental than the observed popu- 
lation itself, it is vital to understand what biases, 
be they observational or astrophysical, may af- 
fect the inferred properties and how they manifest 
themselves. This point is especially consequen- 
tial for high-z Lya-selected studies where typical 
detections are faint and galaxies often go unde- 
tected in the continuum. Lya is a resonant line 
and its formation is strongly affected by a com- 
plex radiative transport. On a physical level the 
regulation and transport of Lya is known to be af- 
fected by dust, the topology and ionization of the 
ISM, Hi kinematics, and viewing geometry. Such 
insights have been gleaned empirically from spec- 
troscopic obsj3ryations^^ 

galaxies dGiavalisco et"al]|l996t[K"untn et aljl998l 
iMas-Hesse et al.ll2003l) , flanked by a sophisticated 
theoretical and computational attack on the prob- 
lem of Lya escape physi cs (iNeufeldl Il990i: lAhnl 
2004 lHansen fc Ohl 120061 IVerhamme et alB oOG: 



Laursen &; Sommer-Larsenl l2007) . 



Since Lya is a resonance line and can be ex- 
pected to be substantially spatially decoupled 
from UV continuum radiation and other nebular 
lines, the picture yielded by UV-targeted spec- 
troscopy is limited. Thus the imaging approach 
becomes an invaluable complement to the previous 
IUE and HST spectroscopic studies. Lya imaging 
at z w is technically possible with HST using the 
Wide Field and Planetary Camera (WFPC2) al- 
though the instrumental throughput at Lya would 
make it very inefficient. A much more economical 
approach would be to use the UV-optimized chan- 
nels of the Space Telescope Imaging Spectrograph 
(STIS) or Advanced Camera for Surveys (ACS) 
which both offer higher system efficiency and bet- 
ter defined and non red-leaking filters. In HST 
cycle 11 we began an imaging study to examine a 
handful of local star-forming galaxies in Lya using 
H ST/ ACS. Firs t resul ts from this were presented 
in iKunth et all (|2003l) . although technical uncer- 
tainties about how to subtract the continuum pre - 



vented a deep analysis and Haves et al. ( 2005f) 



demonstrated the need for additional off-line ob- 
servations to aid in the continuum subtraction. 
This additional data has since been obtained for 
the remainder of the sample and the entire data- 
set has now been processed. Lya line-only images 
are being released to the community (Ostlin et al. 
2009). However, since the process of subtracting 
the continuum in this study is far from trivial and 
future imaging studies depend upon the method- 
ology, this process was deemed to be worthy of 
an article in its own right. We here perform a se- 
ries of tests of synthetic Lya imaging observations 
of low-z targets using configurations available on 
HST and compare different approaches to sub- 
tracting the continuum. 

The article is organized as follows: in Sect.[5]we 
explain the complications involved with Lya ob- 
servations and why conventional continuum sub- 
traction techniques are not applicable to Lya with 
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the current instrumentation; in Sect. [3] we de- 
scribe the methodology and theoretical tests; in 
Sect. S] we present the results and discuss their 
significance; in Sect. [5] we present simulations for 
a study targeting slightly higher redshift and dis- 
cuss some possible augmentation to the method; 
and in Sect. [6] we make some concluding remarks. 



2. Local Lyman-alpha imaging 

2.1. Observing strategies 

The Solar Blind Channel (SBC) of EST /ACS 
offers the F122M filter through which to observe 
restframe Lya at < z < 0.035. This filter peaks 
at 1216A and has a rectangular width (W T ) of 
128. 4A although the pivotal wavelength (A p ) at 
1273. 7A is rather offset from the peak due to the 
red wing. Furthermore, with W r /\ p ~ 0.1, this 
filter cannot be considered narrow compared to 
optical narrowband filters. A number of long- 
pass filters are available through which to sample 
the continuum, the most appropriate for z w 
being F140LP (A p =1527A; W r =252.95A). This 
F122M+F140LP configuration is that used by our 
imaging campaign and is denoted configuration 1 
in this article. Due to the very sharp cut-on on the 
blue side, and the red side being defined by the 
degrading sensitivity of the detector with wave- 
length, all the long-pass filters exhibit near iden- 
tical shapes in the red wing. Adjacent long-pass 
filter combinations could therefore be appropriate 
for Lya imaging of slightly more distant targets 
(e.g. F125LP+F140LP for 0.028 < z < 0.1) and 
this combination is also examined here, denoted 
as configuration 2. Fig. [T] shows the bandpasses 
for both configurations. 

2.2. Continuum subtraction methods 

The art of continuum subtraction hinges en- 
tirely upon scaling of the stellar continuum sam- 
pled by an off-line filter to the bandpass of the 
filter that isolates the line. In the forthcoming dis- 
cussion we make use of a quantity known as the 
Continuum Throu ghput Norma l ization ( CTN) fac- 
tor, first defined in lHaves et all (|2005l) . This is the 
dimensionless quantity that scales the raw count- 
rate in a continuum filter to that expected in the 
on-line filter, accounting for the filter transmission 
profiles, the instrument sensitivity, and the shape 
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Fig. 1. — ACS/ SBC filter combinations for Lya 
imaging. Solid lines represent the filter that trans- 
mits Lya and dashed lines that for the contin- 
uum. Vertical dot-dashed lines arc at A = 1216A, 
the wavelength of both Milky Way Lya absorp- 
tion and the geocoronal emission line. The red 
line shows a synthetic spectrum with artificially 
added Lya. For illustration, a Milky Way absorp- 
tion feature has also been added to the templates 
(log(n H i) = 20.5 cm -2 ). Upper. Configuration 1: 
F122M+F140LP for z < 0.03. Suitable for z « 
although both geocoronal emission and Galactic 
absorption fall in the on-line bandpass and must 
be accounted for. The spectrum is redshifted to 
0.01. Lower: Configuration 2: F125LP+F140LP 
for 0.03 < z < 0.1. Geocoronal Lya is not trans- 
mitted and on-line system throughput is greatly 
improved over the upper example. The spectrum 
is redshifted to 0.06. 
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of the continuum. That is, 



Lya = online — CTN x offline 



(1) 



CTN can be computed with varying degrees of 
complexity. For known or assumed continuum 
slopes it can be computed simply from the in- 
verse sensitivity bandpass characteristics (PHOT- 
FLAM) of the on- and off-line filters. For any 
continuum spectrum, CTN can be computed 
by convolving the spectrum with the instrument 
throughput profiles of the filters, integrating to 
estimate the count-rate, and calculating the ratio. 
This permits features such as spectral shape and 
absorption or emission lines to be accounted for 
in the process. In essence CTN is simply a color 
with a non-standard normalization derived from 
instrument sensitivities. 

There are a number of possible issues that may 
impede the continuum subtraction of Lya. Firstly, 
possible strong stellar Lya absorption may can- 
cel some or all of the nebular emission. This is 
perhaps unlikely to be a strong affect where the 
starburst is very young and dominated by O-stars 
where deep absorption features have yet to de- 
velop. The effect may become very significant 
when Lya is scattered and may be superimposed 
upon older stellar populations with lower effective 
temperature. Indeed, absorption of stellar contin- 
uum by B-stars was initially proposed as an expla- 
nation for the weakness or absence of Lya emis- 
sion in nearby starbursts and for th e early failure 
to de tect Lya emitters at high-z ( Valls-Gabaudl 
1993). Secondly, the Lya feature may be P Cygni 
with some or all of the emission cancelled by ab- 
sorption, especially directly in front of the bright- 
est clusters. Spectroscopic measurements, if of 
sufficient resolution, can isolate the emission seg- 
ment of the P Cygni profile and obtain the emit- 
ted flux, although this is not possible with imag- 
ing. If the absorption segment falls within the 
on-line filter as it will in our case, the local min- 
ima in Lya flux from P Cygni absorption on top 
of clusters can never be corrected for. Notably 
however, this caveat also applies to imaging sur- 
veys at high-z. Ultimately what is called emission 
is a matter of definition; whether it be defined 
as the flux in the emission segment of the pro- 
file, or the flux integrated over both the absorp- 
tion and emission segments. Either way, imaging 
can only assess the latter. Finally, and especially 



for configuration 1, photometry will be directly 
affected by Milky Way Lya absorption and geo- 
coronal Lya emission. Both of these effects can 
be relatively easily accounted for since they are 
not expected to vary strongly across the small an- 
gles subtended by the SBC field-of-view. Geo- 
coronal emission can be removed by background 
subtraction, and Milky Way Lya absorption com- 
puted from the Hi column den sity measured along 
the line-of-sight to the targ et ( Dickey fc Lockmanl 
1990; iKalberla et all 12005) and incorporated into 



the calculation of CTN . 

2.2.1. Purely observational methods 

Typically, for optical emission lines in ground- 
based observations, continuum is sampled using a 
filter positioned as close as possible to, but not 
transmitting, the line in question (or other con- 
taminating spectral features). For example, an 
off-line Ha filter 100A redward of the line corre- 
sponds to AA/A ~ 0.015. Were an intrinsically 
power-law continuum with (3 = —2 assumed to 
be flat (P = 0), this would result in an error in 
the line-center continuum estimate of around 39c0. 
Typically, this would be considered 'good enough'. 

The same is not true for the FUV filter set avail- 
able on HST/SBC. Due to the limited choice of 
filters available, for the F122M+F140LP the on- 
and off-line combination, the pivot wavelengths 
are separated by AA/A ~ 0.2 which would result 
in errors of > 30% if the same misguided assump- 
tion of the continuum slope were made. Errors 
in line-center continuum estimation of this magni- 
tude could translate into severe errors in contin- 
uum subtraction of Lya since the on-line bandpass 
is so broad: if the line does not dominate, then the 
continuum subtracted pixel could quite feasibly 
get the wrong sign (emission could be seen as ab- 
sorption and vice versa) . For a dust-free starburst 
the slope of the FUV continuum is largely con- 
stant (3 ~ —2.6 over the first few 10s Myr, but flat- 
tens quite rapidly (/3 increas es) thereafter, particu - 
larly at higher metallicities ( Leitherer et a,l.lfl999h . 
Since the UV slope is also a strong function of the 
dust reddening (Eb-v), it is not possible to reli- 
ably predict (3 globally, let alone how it varies on 
~ 10 pc scales (i.e. pixel-to-pixel, at SBC sam- 



1 Continuum is assumed to take the form of a power-law in 
A of the form fx oc A* 3 
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pling). One possible alternative would be to sam- 
ple the blue and red sides of the on-line filter but 
there are no filters available for such an observa- 
tion and spectroscopic observations with the STIS 
have shown that the continuum on the blue side 
of Lya is frequently unpredictable, contaminated 
by internal or Galactic absorption features. Flux 
at Lya due to continuum processes must be esti- 
mated from observations on the red-side of Lya 
only. 

The next step would be to take an addi- 
tional off-line observation redwards of the off- 
line filter (e.g. ACS/SBC/F150LP or F165LP, 



( Haves et al. 2005) shows a large diffuse emission 
component which dominates the total Lya lumi- 
nosity. High resolution imaging frequently reveals 
that starbursts are composed of numerous com- 
pact young knots and clusters that dominate the 
luminosity but do no t account for a large fraction 
of the surface area ( Meurer et al. 19951 ). partic- 
ularly at UV wavelengths. It is therefore rather 
likely that the continuum that needs to be sub- 
tracted will not be representative of a region young 
enough to produce Lya, which can emerge super- 
imposed upon regions that appear too d usty for 
Lya to be transmitted (jHaves et alJl2007l ). 



WFPC2/F218W ACS/HRC/F220W, or WFC3/UVIS/F225Vfygure 1 demonstrates inadequacy of assuming 



and extrapolate a power-law continuum to the on- 
line filter. This way CTN would be computed for 
a given spectral slope by convolving the through- 
put profiles with the measured power-law, result- 
ing in a different CTN in each pixel. However, 
this method was demonstrated to result in sig- 
nificant errors in the continuum flux estimation 
when re latively moderate amounts of dust are 
present ([Haves et all 120051 lHaves fc Ostlinl l2006h 
since typical dust reddening modifies the FUV 
continuum in such a way that it becomes incon- 
sistent with the power-law approximation. Fur- 
thermore, such a method provides no estimate of 
stellar Lya absorption which is significant in all 
but very hottest stars. 

Naively, these points could be argued away: 
only the very youngest star-forming regions pro- 
duce enou gh ionizing photons t o generate signifi- 
cant Lya (jCharlot fc Falllll993h over which times 
(3 is essentially constant. High resolution imag- 
ing studies have however demonstrated that in the 
ISM of active starbursts, particularly surround- 
ing massive young clusters, the nebular emission 
can be clearly displaced from the ionizing sources 
(e.g. iMaiz-Apellaniz et al] 12004 ). This is most 
likely due to stellar winds and supernova feed- 
back clearing bubbles in the ISM and is clearly 
visible in some cases in our sample, well exempli- 
fied by ESO338-IG04 (Ostlin et al. 2009). The 
case is further complicated for Lya by its reso- 
nant nature. As a result of multiple scatterings 
in Hi, Lya photons diffuse away from their pro- 
duction sites and are likely emitted from a site 
of last scattering that is not coincid ent with the 



a one-to-one relationship between CTN and (3. 
Here the two quantities are computed for various 
ages and Eb-v by c onvol v ing Starburst99 spec- 



tra ( Leitherer et al. 19991 : Vazquez fc Leitherer 



nebu l ae wh ere they were produced (jHaves et al 
2005L l2007t Ostlin et al. 2009). ESO338-IG04 



20051 ) with the filter t hroughput profiles. For 
this demonstration the ICalzetti et al. I (Il994h at- 
tenuation law is adopted which is frequently used 
to represent dust extinction in starbursts and 
LBGs. The first example shows our own observ- 
ing strategy using SBC/F122M for Lya on-line, 
SBC/F140LP for off-line, measuring j3 between 
F140LP and HRC/F220W. The second example 
shows how, by adopting the F125LP, F140LP and 
F150LP configuration, the degeneracy in CTN vs. 
[3 is significantly reduced. However, with uncer- 
tainties still greater than 40% at certain (3, it is 
still not sufficient for the broad on-line filter. In 
general the extinction laws of the Milky Way and 
Ma gellanic clouds are fo und to be steeper than 
the lCalzetti et al.l (|1994l ) law at given Eb-v over 
the wavelength domain covered by our filters and 
the effect is found to be more severe. 

2.2.2. Methods beyond pure observation 

CTN is the only important quantity to know 
for the continuum subtraction but unfortunately, 
as demonstrated in Section 12.2. 1[ the quantity 
cannot easily be estimated without some knowl- 
edge of the continuum. A solution is compute 
CTN from spectral synthesis models using sup- 
plementary data to constrain the model. Thus 
one must have sufficient information at hand to 
build the right stellar population; stellar age and 
star-formation history, Eb-v, metallicity, and the 
initial mass function (IMF) may all have a signif- 
icant effect. 
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Fig. 2. — CTN vs. (3 for various ages (lMyr 
to 1 G yr) and Eb-v i n the range - 1 for the 
law of ICalzetti et al. ( 1994 ). Lines link points 
of the same age (i.e. follow tracks of dust red- 
dening). Upper. The configuration used for 
our current HST campaign. CTN is generated 
for F122M/F140LP with f3 measured between 
F140LP and F220W. CTN is clearly not a mono- 
tonic function of (3 and degenerate values of CTN 
can be derived for any measured f3, especially for 
bluer colors (e.g. f3 ~ —1). Lower. Same as upper 
but for the configuration 2. CTN is generated 
for F125LP/F140LP with f3 measured between 
F140LP and F150LP. The degeneracy still exists 
although has been quite significantly improved. 



In Haves et al.l ( 2005 ) we demonstrated that 
we could reliably continuum subtract Lya, pixel- 
by-pixel, in our HST observations (configuration 
1) through multi-color spectral modeling using a 
number of additional broadband HST observa- 
tions. In that study four line- free observations 
were used, together with the Starburst99 spec- 
tral evolutionary models. In a,ge-Es-v space we 
showed that we could find non-degenerate solu- 
tions for CTN by sampling (3 (sensitive to both 
age and dust) and the Balmer/4000A break (sen- 
sitive primarily to the age of the stellar popula- 
tion). Grids of multiple color indices and CTN 
were computed from the models for a range of age 
and Eb-v, enabling us to look up the most appro- 
priate value of CTN in each pixel. We also demon- 
strated that several other parameters (metallicity, 
IMF, etc.) had only a small impact upon the reli- 
ability of the method. If x 2 is defined as 



X 



lC- mi (t,EB-v)~d t } 2 -Wi (2) 



where rrii(t, Eb-v) represents the ith model data 
point on the SED for given age (t) and Eb-v, C 
the model normalization factor, di flux measured 
in the ith filter, and Wi the weight for the corre- 
sponding filter, typically the inverse variance. The 
value of C that provides the best fit is determined 
from 

C = -^- 2 (3) 

i 

where C and \ 2 are computed for each possible 
model spectrum, selecting the best models from 
the minimum value of x 2 ■ Figure [3] shows an ex- 
ample x 2 map in t-Es-v space. 

Naturally, if CTN is sensitive to the Balmer/4000A 
break, then a number of contaminants may affect 
the CTN determination, depending upon how 
well the stellar and nebular emission regions are 
resolved. At the earliest times when the nebu- 
lar emission is strong, the Balmer edge may con- 
tribute by bluening the restframe U — B color. On 
the other hand, at later times the 4000A break 
arises from metal line blanketing in late type stars, 
reddening U — B. Thus resolved nebular gas and 
old stellar populations may both affect the SED in 
the optical but not the FUV; not affecting CTN 
itself but its determination should optical SED 



6 



2.0 




Wavelength [ A ] 



Fig. 3. — Left: log(x 2 ) map for a one component stellar fit of age and Eb-v- A degeneracy is clearly seen, 
running diagonally and negatively in age- Eb-v space although note the logarithmic scaling. Data for this 
plot is taken from a diffuse Lya emission region in the vicinity of Knot B of Haro 11. Upper right: Example 



of a multi-component SED fitting to real data from Haro 11 (H aves et al.l 120071 ). Blue shows the nebular 



gas continuum spectrum, the level of which is determined from the Ha observation (note the data point 
at Ha is lies high above the fit, since this observation contains both line and continuum), green shows the 
underlying field-star population of age around 1 Gyr that dominates in the /—band, red shows the current 
starburst with an age of 1 Myr. Black shows the total combined spectrum. Dashed black line shows a 
single component spectral fit. Lower right: Same as upper but extracted from a well resolved Ha shell in 
ESO338-IG04 (Ostlin et al. in prep.). The extremely high Ha equivalent width is apparent since the stellar 
continuum is rather weak. The nebular continuum level is higher than that of the stars in the U and / bands. 
The data-points show the U — B color is negative and the Balmer jump is significant in the total spectrum. 
The nebular component still does not contribute near Lya but demonstrates why control over the nebular 
continuum spectrum is necessary. 
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data-points be relied upon. Since the relative con- 
tribution from stellar populations and nebular gas 
is unknown, and in high-resolution imaging the 
regions may be clearly resolved, it may become 
necessary to account for all these populations in 
our estimation of the true spectrum, and therefore 
CTN. 

Control over the nebular gas contribution may 
be obtained through additional observations that 
directly trace the ionized gas. For given tempera- 
ture in an optically thin nebula, the emission co- 
efficients {a e ** ) may be computed for any allowed 
recombination line, and continuous emission coef- 
ficient (7,,) computed for given wavelength. Thus 
an appropriately strong and well-behaved line (e.g. 
Ha) can be used to estimate the nebular contin- 
uum contribution in each filter, and subtract it 
from each data-point for a given Eb-v ( a g a i n i 
the filters have been selected in order to avoid 
the strongest nebular emission lines). CTN may 
then be computed using the SED-fitting method 
to obtain the stellar-only SED (minimizing % 2 us- 
ing Equation [2]) and reconstruction the composite 
spectrum. 

The current starburst may be superimposed 
upon any underlying stellar population. Since this 
is likely to show a significant 4000A break but lit- 
tle FUV contribution, it may affect the reliabil- 
ity of CTN determination that relies upon age 
fitting. Thus it may be necessary to treat more 
than one stellar population per pixel. Old stel- 
lar populations may dominate the integrated light 
at red wavelengths and such an observation may 
facilitate the decomposition of starburst and un- 
derlying stellar components. Equation [2] can be 
modified to include two populations, a and 6, each 
with differing normalization as 



where 



X 



[C a 



+ Cb ■ m b j 



Wi (4) 



where m a .i and m^j represent the model SED 
data-points, with C a and Cb the normalization fac- 
tors for each population. For a given set of models 
m a and m&, an analytic solution exists for the val- 
ues of C a and Cb that provide the best-fit, as: 



C a 
C b 



S5S2 — S3S4 
- S1S4 

C a S\ 



si 



Si = 

s 2 = 

S3 = 

s 4 = 

s 5 = 



m a .imb,iWi 

i 

22m ati diWi 

i 
i 

^ m b ^diWi 



(6) 



S 2 



(5) 



With the contributions from the two stellar 
populations measured or fit, the total SED can 
be reconstructed by summing the SED in each 
component. Thus CTN is determined from the 
reconstructed model spectrum. 

The fact that stellar ages may be estimated 
provides the first-order solution to another po- 
tential problem: that of the underlying stellar 
Lya absorption feature. For O stars, and de- 
pending on the wind properties, Lya can be ex- 
pected in absorption, filled in P Cygni, or emission 
(net equivalent width s of a few A are expected, 
Klein &; Castor! Il978l ). However, due to the deep 
Galactic Hi column to even the nearest O stars, 
the intrinsic Lya profile is completely absorbed 
and it is likely these predictions can never be em- 
pirically tested. Some B stars however are near 
enough to be observed at Lya and NLTE atmo- 
sphere models appear to be able to match the ob- 
served profile (Lya absorption equivalent widths 
of a few tens of A). Since no observational tests 
are available, it is not deemed meaningful to pur- 
sue this issue further, although it is worth pointing 
out that the stellar feature is accounted for, within 
the limits of current understanding. 

In the following Sections we assess the relative 
power of the various methods described above in 
recovering Lya fluxes and equivalent widths. 

3. Methodology and Tests 

The aim of this study it to test how reliably 
we can recover Lya observables from our imaging 
observations. Ideally we would like to test the vari- 
ous methods of continuum subtraction against real 
observational data for which the true Lya fluxes 
and equivalent widths are known (i.e. spectra). 



S 



Unfortunately there are no real object spectra 
that span the wavelength range between Lya and 
9000A with sufficient spectral resolution, taken 
in consistent apertures for us to test the various 
methods against real data, and we are forced to 
rely only upon synthetic input sources. 

In this study we present a number of com- 
putational tests to assess the performance of 
various methods. To this end we generate a 
set of template spectra usi ng various combina- 
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tions of Starburst99 models ([Leitherer et al 
Vazquez fc Leithererl 12005t h for various stellar 
populations (or combinations thereof), Eb-v, an d 
modify them by adding Lya lines of chosen equiva- 
lent width. We then convolve the template spectra 
with the HST bandpasses to generate synthetic 
SED data-points and feed these into our SED- 
fitting and continuum subtraction software. Thus 
we have full knowledge of the intrinsic spectrum, 
and we are able to compare our output results 
directly to the known input values (CTN, F^ ya , 
Wi,y a ), and test the ability of our fitting meth- 
ods to recover ages and reddening in the stellar 
population(s). 

We simulate noise by randomizing each SED 
point for a given error and, adopting a Monte- 
Carlo approach, we compute various statistics of 
the recovered distributions. For the real HST data 
we employ adaptive binning techniques to obtain 
a minimum threshold S/N, but the choice of this 
value requires testing. Therefore by performing 
these simulations we are able to test the results 
obtained against input S/N, and directly test the 
optimum threshold. 

The method of generating the 'real' input spec- 
tra is outlined in Section 13.11 The various meth- 
ods of subtracting the continuum are described in 
Section 13. 2\ and the actual quantitative tests are 
described in Section l3~3l 



3.1. Input SED generation 

We have previously discussed the fact that SED 
is likely to be the sum of the contributions from 
a current episode of star formation (starburst; 
sb), an underlying component of field stars (fs), 
and emission from nebular gas (neb). All tem- 
plate input spectra are generated from a combina- 
tion of these three components, scaled to a given 
normalization relative to the starburst at 4500A 



(~ B— band); these normalizations are assigned 
the ri{ s and n ne b for the field stars and nebular 
components, respectively. The templates thereby 
consist of two stellar components of variable ages, 
and a single gas component, all contributing differ- 
ent fractions of the B— band luminosity. For sim- 
plicity, the two stellar components are generated 
from the same metallicity and initial mass func- 
tion. The Ha luminosity is computed from the 
nebular gas continuum flux density and applied to 
the gas spectrum before addition. To complete the 
restframe spectrum, it is reddened using a given 
extinction law and Eb-v ■ Finally the spectrum is 
redshifted. Further information on the generation 
of input spectra can be found in Section 13.31 and 
Table [U 

For each template spectrum, we generate a 
set of 'observed' SED data-points by convolv- 
ing the spectra with the HST filter profiles. 
For configuration 1 (our current imaging cam- 
paign) the complete filter list is: SBC/F122M 
(Lya on-line), SBC/F140LP (~ 1500A con- 
tinuum), HRC/F220W (~ 2200A continuum), 
HRC/F330W (~ ?7-band), WFC/F435W (~ 
B-band), WFC/F550M (medium V-b&nd, serv- 
ing as line- free continuum filter near Ha), WFC/FR656N 
(linear ramp narrow-band filter centered upon 
restframe Ha), and HRC/F814 W (~ /-band). 
For configuration 2 the filter set replaces F122M 
with F125LP while F140LP remains the FUV 
continuum filter as shown in Figure [1] For the 
optical component of this configuration we sample 
the UV/optical continuum using WFPC2/F336 W, 
F439W, and F814W, using the fixed narrow-band 
filter F673N to observe redshifted Ha. 

3.2. Continuum subtraction methods 

For each SED we apply various methods of es- 
timating and subtracting the continuum and ex- 
amine how well they return the known restframe 
quantities that were input. We identify five possi- 
ble methods of estimating CTN and subtracting 
the continuum at Lya: 

I Assuming the slope of the continuum between 
the off-line and on-line filters to be flat in f\. 
I.e. CTN is the ratio of the HST PHOT- 
FLAM values for each filter. 

II Assuming the continuum takes the form of 
a power-law in fx (oc A' 3 ) and extrapolating 
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the slope as measured between the NUV fil- 
ter at ~ 2200 A and the off-line FUV filter 
(F140LP). 

III Fitting a single Starburst99 spectrum (com- 
bined stellar+nebular components) to the 
continuum SED points (i.e. excluding the 
Ha observation) using Equations [2] and [3] 
Then computing CTN from the best-fitting 
(age and Eb-v) spectrum. 

IV Continuum-subtracting Ha and using the Ha 
flux to reconstruct the SED due to nebular 
continuum processes. Then fitting a single 
stellar-only Starburst99 spectrum (again fit- 
ting age and Eb-v) to the real SED data- 
points with the nebular contribution sub- 
tracted (i.e. real data-points - nebular data- 
points) using the same method as method 
in. 

V Same as method iv but fitting 2 stellar com- 
ponents to the nebular-subtracted SED using 
Equations 14151 and[U In this case, the fitting 
is performed over all ages with age allowed to 
vary in both stellar populations. 

3.3. The tests 

The most important returned values are, of 
course, -FLya and W-Lya, although from methods 
where SED fitting is employed, age(s), normaliza- 
tion factor(s), and Eb-v are also returned. Mon- 
itoring these returned values in addition to those 
relating to Lya permits a deeper examination of 
the performance of the SED fitting. 

We begin by defining a fiducial template spec- 
trum, denoted fid. This is generated from what is 
thought to be a typical starburst, capable of pro- 
ducing Lya with moderate extinction and metal- 
licity and a Salpeter IMF. For the starburst, a sin- 
gle stellar population of age of 5 Myr is selected 
since ionizing photons are a requirement for the 
production of Lya. According to the Starburst99 
template spectra, the nebular gas contribution at 
4500 A (n nc b, defined relative to the stellar spec- 
trum of and unresolved point-source) is ~ 0.05 
at this age, and this value provides the fiducial 
n nc b- For the fiducial model we add a moder- 
ate field-star population with age (arbitrarily se- 
lected) 5 Gyr, scaled to give a contribution of 
50% the starburst luminosity at 4500A (rif s =0.5). 



The fiducial reddening for the composite spectrum 
was selected to be Eb-v= 0.2 using the SMC 
law. The value of Eb-v is an approximate mid- 
point in the measured extinct ion for most of our 
galaxies fsee lAtek et al. I l2008l) . The choice of law 
is motivated by the fact that the SMC law ap- 
pears to provide the best fits to the resolved spec- 
tra of blue compact and irregular galaxies (e.g. 



Mas-Hesse fc Kunthl 11999). Since Wj, ya can es- 



sentially take any value between damped absorp- 
tion and super-recombination values, we add Lya 
lines with a range of equivalent widths: -50A to 
represent the deep absorption seen in some lo- 
cal objects (e.g. iZwl8); OA, corresponding to 
no Lya feature; 10A to approximate the global 
values measured in s ome of low-z objects (e.g. 
Giavalisco et al.|[l996h : and 100 A corresponding to 



high-z Lya-bright galaxies or diffuse emission re- 
gions in resolved objects at z r* 0. We define a 
number of modifications to all of the ingredient 
parameters in the input spectrum which are listed 
in Table [1] 

Firstly, all tests are performed at infinite S/N 
to test the reliability of the code by insuring all 
the input parameters are returned. 

Observational S/N per resolution clement 
varies significantly with position as a result of 
natural morphological variation in surface bright- 
ness, and band-to-band as morphology differs with 
wavelength. In our ACS data-set we make use of 
adaptive binning to bin together pixels until a 
threshold S/N has been met, conserving surface 
brightness in each conglomerate bin ( "spaxel" ) . 
Without investigation it is not possible to know 
what S/N is required in order for the various 
techniques to become optimal; either converging 
on the input values or to some systematically off- 
set values. By varying S/N in our simulations, we 
are able to test this directly. To each data-point 
in the input SED (see Section 13. ip , we assign an 
error, based upon the S/N we want to test. Us- 
ing this error we then regenerate 1000 SEDs using 
Gaussian-variates and apply all our continuum 
subtraction methodologies to each one. Perform- 
ing such Monte-Carlo simulations with increasing 
S/N allows us to investigate the optimum S/N 
required to return mean values consistent with 
the those input, or find convergence of the mean 
to systematically offset observables. For each of 
the model galaxy SEDs defined in Table [1] we run 
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Table 1 

Parameters for SED generation 



Parameter 


Unit 


Fiducial value 


Non-fiducial / Range 




Code 




A 


10 


-50, 0, and 100 


t flcLwlO; fid.w-50; fid.wO; fid.wlOO ] 


Eb-v 


mag 


0.2 


0.0 and 0.5 


[ ebvCLw . . . ; 


ebv05_w. . .] 


Reddening law 




SMC 


Calzetti et al. f 1994) 


[ cal_w . . . ] 




Starburst age 


Myr 


5 


20 and 100 


[ agesb20_w. 


. ; agesbl00_w. . .] 


Field star age 


Myr 


5000 


200 


[ agefs200_w 


. . ] 






0.5 


and 5 


[ nfsO.w. . . ; 


nf s5_w. . .] 


"nob 2 




0.05 


and 1 


[ nneb0_w . . . 


nnebl_w . . . ] 


IMF a 




-2.35 


-1.85 and -2.85 


[ imf 185_w. . 


; imf285_w. . .] 


Mctallicity 


Z 


0.08 


0.001 and 0.040 


[ met sub _w. . 


; metsuper_w . . . ] 


Stellar atmosphere 3 




Atm5 


Atm3 and Atm4 


[ atm3_w . . . ; 


atm4_w . . . ] 


Redshift 




0.01 









Note. — Model parameters in the third column constitute the fid_wl0 fiducial template spectrum case, with the 
non-fiducial values for W-^ ya constituting standard cases of f id_w— 50, f id_w0, and f id_wl00. 

1 defined relative to Starburst population at 4500A (i.e. ~ relative B— band luminosity). 

2 defined relative to Starburst population at 4500A. 0.05 corresponds to typical default unresolved nebular fraction at 
5Myr from Starburst99. and 1 correspond therefore to zero nebular contribution and a 'boost' by a factor of ~ 20. 

3 Designated Starburst99 codes. Atm3 is Lejeune atmospheres for stars with plane-parallel atmospheres and Schmutz 
atmospheres for stars with strong winds. Atm4 is like Atm3 but with Hillicr atmospheres for stars with strong winds. 
Atm5 is like Atm4 with Pauldrach models for O stars. 



1000 Monte-Carlo iterations at each S/N with 
S/N varying between 1 and 50. Important statis- 
tics of -FLya and W^ ya ; the mean and median, 
standard deviation, and skewness of the derived 
distribution are retained. 

Firstly we set S/N to be the same in all band- 
passes; if the worst quality band is chosen to gen- 
erate the binning pattern, S/N in the correspond- 
ing bins in other bands should have S/N exceed- 
ing the threshold. However, due to differing mor- 
phologies, this is unlikely to be the case in gen- 
eral, and some bins, in some images, may fall short 
of the desired S/N threshold. We assess the im- 
pact of this to the overall fitting by systematically 
dropping S/N to zero in individual bands for the 
fiducial SED. This allows us to examine the ro- 
bustness of returned Lya-related quantities in the 
case where lower-quality data have been obtained 
in certain bandpasses. Similar tests are then per- 
formed dropping S/N to zero in adjacent pairs of 
filters simultaneously. 

To test the reliability of the methods for a wide 
array of model galaxies we augment the fiducial 
SED, one parameter at a time, always running the 
SED-fitting and continuum subtraction code us- 
ing the standard Starburst99 parameters and SMC 
reddening law. Age and reddening, of course, are 
always the free parameters in the fitting (both age 



components are fit for method v). This way we 
can assess the reliability of our results at differ- 
ent stellar ages, Eb-v, an( A with different con- 
tributions from the underlying stellar population 
and nebular gas. The same is also true when we 
test metallicity, IMF, reddening law, enabling us 
to test the impact of failing to select the correct 
values for the fitting. The effect of metallicity is 
found to be have sufficient impact that additional 
tests are performed, see Section H.2.31 The param- 
eter space covered is listed in Table [1] Parameter 
dependency tests are performed for both observa- 
tional configurations with a detailed discussion of 
the results for configuration 1 , including the effects 
of poor S/N in certain filters, presented in SectHJ 
A summary of the results for configuration 2 is 
presented in Sect El 

4. Results and discussion 
4.1. The fiducial model 

In this subsection we first present and discuss 
the results obtained for the fiducial SED with the 
same S/N in all bandpasses in Sect 14.1. II followed 
by results where single and adjacent pairs of filters 
have S/N = in Sect l4~l~2l 
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Fig. 4. — Statistics of the returned distribution from 1000 Monte-Carlo iterations on the fiducial setup with 
Whya of 10A. Top panels show the mean, center the la standard deviation, and lower the skewness. Left: 
Continuum subtracted Lya flux in normalized units (i.e. returned mean divided by input flux). Standard 
deviation is also given in normalized units. Right: Lya equivalent width. Colors, line styles, and weights 
show the various methods of estimating the continuum and arc labeled in the legend. 



4- 1-1. Equal S/N in all bands 

Figure 2] shows the returned statistics for the 
fiducial model as a function of S/N . Statistics 
shown include the mean (upper), standard devia- 
tion (center), and skewness (lower). Firstly, these 
plots demonstrate that simple assumptions about 
the continuum slope (either f3 = (method i) 
or (3 extrapolation (method n)) both result in an 
overestimate of CTN . Therefore the continuum is 
over- subtracted and a 10A emission line is seen as 
an absorption feature with WL ya ~ — 40A. This is 
the result of: (a) the large offset in A between the 
off- and on-line filters; (b) the fact that the ~ 100A 
wide on-line bandpass is far from line-dominated; 
and (c) the modest value of Eb-v~ 0.2 is suffi- 
cient to reduce the total flux in F122M below that 
of F140LP. It is clear that a good understanding 
of the continuum is essential, not only between 
the filters but across the on-line bandpass itself, 
and that spectral modeling of some level of so- 
phistication is a requirement. It should be noted 
that at S/N — 10 in the individual bandpasses, 
S/N in the continuum subtracted Lya distribu- 
tion is only around 0.5, although obviously this 



improves with increasing W^ya as the line starts 
to dominate. For the ficLwlOO case, S/N — 4 
is seen in the returned Lya flux distribution for 
S/N = 10 in all filters. Of the three SED-fitting 
methods, technique v slightly out-performs ill and 
iv, thanks to its inclusion of treatment of the un- 
derlying stellar population and nebular gas, even 
though they are only minor contributors to the 
FUV flux for the fid SED. It should also be noted 
that all the continuum subtraction methods show 
a positive skew in the Wi^ ya distribution, even at 
S/N > 30. This is due to fact that W is a ra- 
tio, and that the inverse of a Gaussian distribu- 
tion always shows positive skewness - this should 
be present in all equivalent width estimates, irre- 
spective of data-set or observational methodology 
( for a discussion seelDaw s on et al' . 2004; iHu et al 
20041 : lHaves fc 6stlinll200d ). 



An estimate of the age of the stellar popula- 
tion also provides, in part, the solution to another 
potential problem, that of the unknown underly- 
ing stellar absorption at Lya. While still poorly 
tested, Lya features in the models still provide the 
best estimate of the underlying absorption avail- 
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able and are included in computed values of CTN 
and therefore the estimate of the continuum flux. 



4-1-2. Reduced filter sets or poor S/N 

The upper panel of Figure \5\ shows the effect of 
removing a single filter from the fitting routine for 
the f icLwlO SED. Naturally, only the three con- 
tinuum subtraction methods that employ SED fit- 
ting are shown. We now show the mean Wt, ya ob- 
tained with the filter removed, normalized by that 
obtained with all filters included, with the same 
S/N in all included bandpasses. This normalized 
distribution has clearly been shown to converge in 
flux at S/N ~ 5 (Figured]). The results are nois- 
ier than those shown previously for two reasons: 
firstly because the fits are intrinsically noisier due 
to one fewer data-point being present, and sec- 
ondly because of added noise from the normaliza- 
tion. The example shown illustrates the removal 
of the F220W filter, although results are largely 
indistinguishable when other single filters are re- 
moved. This is because we have a maximum of 
5 model parameters to fit and in the case where 
a single filter is removed, we still have sufficient 
data-points remaining to avoid a degeneracy. Pro- 
vided we obtain the minimum S/N threshold of 10 
as shown in Section 14.1.11 in 5 of our 6 continuum 
bandpasses, we can feel safe about the recovery of 
Lya fluxes. 

The situation changes, however, when two fil- 
ters are removed from the fit as can be seen in 
the lower panel of Figure [5] This plot shows the 
same as the upper panel when the two UV filters 
(F220W and F330W) are assigned S/N = 0. The 
two single stellar component fitting methods (in 
and iv) now consistently overestimate Wt, ya by 
around 20% at S/N ~ 10 which actually becomes 
worse as S/N increases in the other filters. This 
is due to more consistent selection of the wrong 
stellar parameters by the code and the example 
has been selected to show poorly recovered observ- 
ables; the loss of both UV filters is the most detri- 
mental. Without either F220W or F330W, we 
have no sampling of the FUV continuum slope or 
the Balmer/4000A break, essentially resulting in 
the recovery of any ages and Eb-v- Clearly main- 
taining S/N ~ 5 sampling of the Balmer break is 
preferable to obtaining extremely high S/N ob- 
servations in other bandpasses; as discussed in 
Wik lind et al. I (|2008l) the Balmer break is instru- 
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Fig. 5. — Average values of W\, ya returned by 
the SED fitting software for the fiducial model 
galaxy with WLya = 10A when datapoints are re- 
moved from the fitting procedure. S/N is the 
same in all the remaining bandpasses. W^ ya is 
now normalized by its value when all filters are 
included, where convergence was previously found 
after S/N ~ 5 (see Figured]). The upper panel 
shows the removal of the F220W filter. The re- 
moval of both F220W and F330W is shown below. 
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mental in resolving the degeneracy between age 
and reddening. We cannot conclude that methods 
in and iv out-perform v in this case but we can be 
quite certain that we need to reach the threshold 
S/N in 5 of our 6 continuum bandpasses in order 
for any of the SED-fitting methods to yield robust 
results. Four filters is not deemed to be sufficient 
for method v. 

4.2. Parameter dependencies 

We now assess the impact of modifying the var- 
ious parameters that go into the construction of a 
the composite galaxy spectra. The fiducial pa- 
rameters are always used internally for the fitting 
(third column in Table Q}, fitting only ages and 
Eb-v- 

4-2.1. Stellar ages and dust reddening 

Table Q] shows the various parameter modifica- 
tions made to the fid galaxy template for star- 
burst age, field star age, and Eb-v ■ For the SED- 
fitting methods (in, iv, and v) results are found 
not to be be significantly discrepant from those 
presented for the fiducial case. In fact, aside from 
details in the noise, the flux and equivalent width 
plots are indistinguishable from those in Figured] 
This is not the case for methods I and n which 
vary wildly with starburst age and Eb-v- Natu- 
rally the behavior of the modeling methods is to 
be expected since age and Eb-v are our primary 
fitting parameters and the effect of dust can be 
well controlled in the event that the chosen ex- 
tinction law well describes the intrinsic deviation 
from the dust-free starburst. 

Exchanging the SMC extincti o n law for the 
attenuation law of Calzetti et al. ( 1994 ) is more 



detrimental: these two curves are not near equiv- 
alent over the wavelength domain we are sam- 
pling. For the cal_wlO case, techniques in and 
iv result in convergent W\, ya estimates of 75 and 
112 A respectively while method v converges at 
WLya=20A by S/N = 10 (still an overestimate 
by a factor of 2). For strong emission (cal_wl00), 
method v overestimates W\, ya by just 10% at 
S/N — 10 while in and iv still consistently fail by 
a factor of 2. However, clearly an accurate model 
of the internal extinction law in the target galaxy 
is a key parameter in the method. 

Several solutions to this exist. Firstly, the red- 



dening law itself could be incorporated as a free 
parameter, only increasing computation time by 
a factor of a few. Including individual extinction 
laws (e.g. SMC, LMC, or Galactic curves) in the 
SED fitting is physically rather poorly motivated 
and a better alternative may be to base the choice 
of law on observation. For example FUV spec- 
troscopy of the objects could indicate the presence 
or absence of a 2 175 A gr aphite feature to moti vate 



the choice of cur ves (see Puget fc Leger 19891 and 



the discussion in Mas- Hesse &: KuntI 



;gerilib)e 
hl ll999h 



4-2.2. Nebular gas and underlying stellar popu- 
lation 

Figure [6] shows the average returned W\, ya as a 
function of S/N when the nebular gas component 
contributes equally with the starburst at 4500A 
(nnebl_wlO; left), and when the underlying pop- 
ulation contributes 5 times that of the starburst 
(nf s5_wl0; right). The left plot demonstrates how 
a single component fit (in) breaks down in the 
nnebl_wlO model when the nebular gas compo- 
nent is dominating the SED (i.e. regions where 
wind-blown Hn shells or filamentary structure are 
resolved). Method in converges with a returned 
W\, ya of around 40A in this case, and a similar 
arithmetic (+30A) overestimate for all the input 
equivalent widths. This demonstrates the need for 
an independent measure of the nebular gas spec- 
trum as methods IV and v well recover Wj^y a . 

The right plot shows the necessity to also con- 
trol the underlying stellar population when fit- 
ting stellar ages. The ordinate axis does not 
show methods in or iv which both converge at 
W^Lya~ 200A. This corresponds to an overestimate 
of an order of magnitude, with overestimates seen 
for all input Why a - Notably, W\, ya converges at 
around 400A in the nfs5_wl00 case and almost 
+100A in the nf s5_w-50 case using methods ill 
and IV. 

4-2.3. Metallicity 

Two additional metallicities were initially 
tested: the minimum sub-solar value (Z = 0.001; 
met sub) and maximum super-solar value (Z = 
0.040; metsuper). The upper panel of Figure [7] 
shows the returned values of W^ya for the mod- 
els with modified metallicities metsub_wl0 and 
metsuper_wlO, again normalized by the returned 
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signal-to-noise signal-to-noise 

Fig. 6. — As in the right panel of Figure 0] but with different input parameters. Left: nnebl_wlO in which 
the starburst and nebular gas component contribute equally at 4500A. Methods I and n both converge 
significantly off the lower end of the ordinate axis. Right: nfs5_wl0 in which the field star population 
outshines the starburst at 4500 A by a factor of 5. Methods III and iv both converge significantly off the 
ordinate axis at equivalent widths of almost 200A. 



values for the f id_wl0 model as described for Fig- 
ure [5] Interestingly while one represents a de- 
crease in metallicity and one an increase, both 
modifications have a similar impact on the result- 
ing observables: the continuum flux in the on-line 
filter is consistently underestimated. The under- 
estimate is around 5% in the both cases but, for 
the weak emission cases this is enough to cause 
overestimates in Whya of 10A (a factor of 2) using 
method v and worse for III and iv. The effect 
is less noticeable in the cases of stronger Lya 
emission where the underestimate of the on-line 
continuum flux linearizes and converges to over- 
estimates proportional to the error on the contin- 
uum. In a similar manner to when filters were 
removed from the fit (Sect. E3~l~2")) . increasing S/N 
means that unrepresentative fits are found more 
consistently, resulting in less spread in the output 
values than at low S/N where the spread is the 
result of noisy fits. In the case of metallicity, poor 
recovery of stellar age is found to be the culprit 
behind deviant Lya observables due to the impact 
of metallicity on the 4000A break. 

The apparent strong effect of unconstrained 
metallicity appears to be a concern and has 
prompted further tests. In the second round 
of testing, all available metallicities were exam- 



ined; Z = [0.001,0.004,0.008,0.020 and 0.040]. 
A set of template spectra was generated for each 
one, and each model galaxy was run through the 
continuum-subtraction software five times, once 
for each metallicity. As demonstrated, quite sig- 
nificantly inaccurate results can be returned if the 
contrast between real and assumed metallicity is 
large (factors greater than around 2). However, 
using method v and for all cases with weak Lya 
emission (WLy Q =10A), the discrepancy in Lya 
fluxes is reduced to less than the returned statis- 
tical error-bar provided that (a) metallicities are 
not discrepant by more than one adjacent step in 
the grid, and (b) S/N = 5 or greater has been 
obtained in all bands. Fortunately very strong 
metallicity gradients are not often found across 
starburst regions and a single long-slit spectrum 
is likely to provide sufficient information. 

4-2.4- The initial mass function 

The lower panel of Figure [7] shows the normal- 
ized returned values of Wi, ya for the models with 
modified IMFs: the flatter imfl85_wl0 and the 
steeper imf 285_wl0. IMFs are modified in both 
the starburst and underlying stellar components 
for the generation of the template spectra. Meth- 
ods in and IV here appear to perform in a simi- 
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Fig. 7. — As in Figure [5] for different metallici- 
ties {upper) and initial mass functions (lower). In- 
put Why a is lOA. SED-fitting methods only are 
shown and methods are labeled. Wr, ya is nor- 
malized by the resulting values obtained from the 
fiducial model, output values of which have previ- 
ously been shown to converge at around S/N ~ 5. 
For metallicity, the thick lines represent sub-solar 
and the thin lines super-solar metallicity. Re- 
garding the IMF, the thick lines show the flatter 
imf 185_wl0 model while the thin lines show the 
steeper imf 285_wl0 model. 



lar manner to when metallicity is modified: nor- 
malized Whya deviates further from unity with in- 
creasing S/N as poor fits are consistently found. 
However, the two-component stellar fit method al- 
ways recovers a mean Wz, ya that scatters around 
the 1-line. This results from the fact that IMF 
and star-formation history both alter the current 
stellar mass distribution: fitting two components 
appears to mimic the effect of a singly modified 
IMF. 

4-2.5. Stellar atmosphere 

Changing the UV stellar atmosphere model 
makes no appreciable difference to any of the 
methodologies. All SED-fitting methods appear 
to return W^ ya consistently around the desired 
value and results are indistinguishable from those 
presented in Figure [4] While testing one stellar 
atmosphere model against another in this manner 
may not be greatly meaningful, it is at least re- 
assuring that the choice of atmosphere does not 
impact upon the recovered values. The stellar at- 
mosphere models differ quite significantly in their 
ionizing output and in the detail of line features. 
However the ionizing output only affects the neb- 
ular fluxes which are treated independently, and 
the discrete line features have only a minor im- 
pact upon integrated colors, even in the FUV. 

5. Future studies and possible further im- 
provements 

So far we have only discussed simulations re- 
lating to the data-sets we have already obtained, 
targeting the lowest redshifts possible. However, 
this redshift severely limits the number of poten- 
tial target galaxies, and any future Lya imaging 
studies will require larger volumes, especially in 
order to target the more luminous analogues of 
high-z star-forming galaxies. Beyond z as 0.03 
the optimal choice of filters changes, and can be 
performed using ACS/SBC and adjacent pairs of 
long-pass filters. Here we present some simulations 
for such a study. We also discuss some method- 
ological alternatives using augmented datasets. 

5.1. A potential local Lya imaging study 
with HST 

As discussed in Sect. 12.11 and illustrated in 
Fig.m we adopt the combination of F125LP/F140LP, 



16 



restricting us to the broad redshift range of 
0.028 — 0.09, although any adjacent pair of long- 
pass filters can be used. With the current inactive 
status of ACS CCD channels and uncertain fu- 
ture of WFC3, we opt for bandpasses available 
on WFPC2 to cover the optical domain, although 
this setup could easily be ported to configura- 
tions on both these cameras. Filters chosen for 
the continuum and Ha observations are listed in 
Section. 13. li lt is also worth noting that with this 
configuration, we have no medium-band line-free 
filter near Ha and the continuum subtraction of 
Ha must rely on interpolation between F439W 
and F814W. Simulations presented here adopt 
the redshift z = 0.029 so as to shift Lya into 
the on-line bandpass without shifting Ha out of 
F673N. 

Figure [H] shows some example results from sim- 
ulations using configuration 2. They are the same 
results as shown for configuration 1 in Figure [HI 
These Figures show the recovered Lya equivalent 
width from the nebular gas dominated and field- 
star dominated templates, and are highly resem- 
blant of those presented for configuration 1. Con- 
figuration 2 has the slight advantage over 1 of us- 
ing two filters with near-identical red wings, which 
isolates a well-defined on-line bandpass shifted 
slightly nearer to the pivotal wavelength of the 
continuum measurement. However, configuration 
2 does still suffer from the same drawback of a 
wide on-line bandpass and, for S/N = 10 in all 
bands, S/N in the continuum subtraction is still 
around 0.5 for the ficLwlO model, the same as 
configuration 1 (Section l4.1.ip . 

Configuration 2 does include two fewer contin- 
uum filters than 1, which was previously shown 
(Section 14.1.21 and Figure [5]) to be inadequate: 
the lower panel shows how Wi,y a is very poorly 
recovered when the NUV and [/—band filters 
are removed. This was designed to exemplify 
the poor recovery of Lya observables and this 
case was especially bad due to the loss of both 
UV filters, and no sampling of the FUV slope 
or Balmer/4000A break. The inclusion of both 
F220W and F 330 W for configuration 1 is likely 
to include some level of redundancy. The contin- 
uum filters of configuration 2 are ideally placed 
and the most important spectral features re- 
main well sampled: F140LP+F336W samples 
f3, F336W+F439W samples the Balmer/4000A 



break, and F439W+F814W simultaneously allows 
for the continuum subtraction of Ha and the con- 
straint of the field star population. In light of 
these results and those presented in Section 14.1.21 
we determine this to be the the absolute minimum 
requirement in SED coverage. 

5.2. Stellar and nebular extinction 

One area in the methodology that could be 
identified as a weakness is the treatment of red- 
dening: we treat stellar (Eb-v,*) and interstel- 
lar extinction (Eb-v.is) m the same, possibly 
sub-optimal, way: by locking them together. 
Eb-v,* and Eb-vjs are known to differ quite 
substantially in some cases, with derived values 
of E r-v.* lower tha n Eb-v,is ([Fanelri et al. 1988; 



Calzetti et al.l fl993 ). This seems to be a geo- 
metric effect due to winds from massive stars ex- 
pelling the ionized ISM and dust, reducing the 
extinction derived from the stellar UV slope and 
concentrating the du st in the Hil shells and fila - 
mentary structures ( Maiz-Apellaniz et al. 1998). 



Mas-Hesse fc Kunthl (fl999h found that the dis- 
crepancy grew with increasing age over timescales 
consistent with stellar evolution and supernova 
enrichment. It is frequently the case that we can 
resolve nebular structure down to the resolution 
limit, although we will can never know the geome- 
try along the line-of-sight. How the dust in front 
of stars and in ionized regions combine depends 
upon the unresolved ISM geometry and it may, 
in such cases, be preferable to treat nebular and 
stellar reddenings independently in the method. 

Decoupling could be achieved simply in the 
fitting procedure, by introducing Eb-v.is as an 
extra fitting dimension that applies only to the 
nebular SED, allowing both values of Eb-v to 
vary independently. However, the number of SED 
data-points is not sufficient for the inclusion of 
an extra degree of freedom, and currently we 
deem this too computationally expensive to in- 
clude in the fitting algorithm when we have > 
10 6 pixels per image. While an empirical re- 
lationship has been present e d betw een the two 



quantities by ICalzetti et al. (2000): Eb 



-v,* 



(0.44 ± 0.033) Eb-v,is, it was derived for a sam- 
ple of galaxies observed with much poorer resolu- 
tion and, in individual HST pixels, the spread be- 
tween the quantities is likely to be so large that the 
two quantities are completely decoupled. Alterna- 
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signal-to-noise signal-to-noise 

Fig. 8. — Example results using observing configuration 2. Input Wh ya is 10A. Left: nebular emission- 
dominated region (nnebl_wlO). Again, methods I and n converge well below the lower limit of the ordinate. 
Right: region dominated by and old stellar population (nf s5_wl0). Methods in and IV converge above the 
upper limit of the ordinate. All continuum subtraction methods are shown and color-coded as in Figure [6] 



tively, Eb-v,is m &y be measured directly by using 
an emission line decrement. Typically the Balmer 
decrement Ha/H/3 would be used where H/3 could 
also be obtained from HST using linear ramp or 
selected narrowband filters. A further possibility 
would be to use a NIR emission line (e.g. Paa or 
Br7) observed from the ground using a large adap- 
tive optics imager. Indeed, an independent evalu- 
ation of the nebular reddening is precisely what is 
required for an appropriate astrophysical compar- 
ison with Lya. A direct measurement of the inter- 
stellar reddening would permit us to lock Eb~v,tb 
on the nebular gas component and fit Eb~v,* to 
the stellar continuum. To this end we are in the 
process of observing a subset of our current sam- 
ple in the Br7 line. Such methods will be tested 
in a forthcoming study. 

Ultimately, however, the continuum we need 
to estimate and subtract is predominantly stellar, 
and it may simply be that treatment of redden- 
ing on the stellar continuum is all that is required. 
Further dedicated tests will be presented when the 
data have been acquired and processed. 

6. Conclusions and Summary 

Using synthetic spectra of starburst galaxies 
we have examined various methods of producing 
continuum subtracted line-only Lya images using 



currently available imaging modes on the Hub- 
ble Space Telescope. We have assessed and com- 
pared various methods of continuum subtraction 
that vary in their complexity and attention paid to 
possible behavior of the continuum. We have pre- 
sented examples covering a wide array of starburst 
parameters for two observational configurations: 
our ACS/SBC imaging campaign of local star- 
bursts (online/offline =F122M / F140LP), and 
for studies using adjacent pairs of long-pass filters 
(F125LP / F140LP), targeting slightly higher-z. 
Our main conclusions are: 

• Making simple assumptions about the shape 
of the far ultraviolet continuum slope (/3; e.g. 
assuming its slope or extrapolating the slope 
from observations on the red side only) leads 
to estimates of the Lya flux and equivalent 
width that are seriously discrepant with the 
true values. Some spectral fitting is shown 
to be essential, even for the most basic cases. 

• In our methodology we fit only the age of 
the stellar population and the dust redden- 
ing, requiring at least datapoints on the SED 
that sample the UV continuum slope and 
Balmer/4000A discontinuity. All other pa- 
rameters (extinction law, metallicity, IMF 
etc.) take assumed values; we then inves- 
tigate the systematic effects incurred when 
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these quantities differ from their true val- 
ues. 

We need to bin together pixels until signal- 
to-noise of between 5 and 10 has been ob- 
tained. 

Age determinations may be contaminated by 
boosted nebular gas emission and an under- 
lying stellar population. These can be ac- 
counted for by constraining the nebular gas 
contribution using an estimate based upon 
the Ha emission flux, and contribution from 
old stars by fitting multiple stellar compo- 
nents, respectively. This requires a data- 
point redwards of Ha, both for the contin- 
uum subtraction of Ha and an estimate of 
the contribution from field stars. We have 
determined the /—band to be functional. 

The initial mass function is shown not to 
be a necessary parameter to include in 
the fitting process and results are largely 
unchanged when non-standard IMFs are 
tested. Metallicity has a much more sig- 
nificant impact upon the recovery of the 
Lya flux, since it changes the rate of stel- 
lar evolution and reddens stellar continua. 
This requires an independent determination 
of the metallicity. The choice of stellar at- 
mosphere model has no discernible impact 
upon the recovered Lya observables. 

If the metallicity and reddening law are 
known quantities and S/N — 5 has been ob- 
tained in all bandpasses, no significant sys- 
tematic effects are seen in our continuum 
subtracted fluxes or equivalent widths. In- 
creasing S/N to 10 is shown to significantly 
reduce the scatter. For input WLy Q = 10A, 
we are able to recover Lya fluxes accurate 
to within 30% of the true value for all tested 
parameter space. This improves to better 
than 10% for stronger Lya emission with 

w LyQ =iooA. 

We have also presented simulations for a 
very similar study that uses adjacent com- 
binations of SBC long-pass filters to isolate 
Lya. Due to the near-identical red wings, 
this could naively be thought to mitigate 
many of the issues surrounding continuum 



subtraction. However, due to the braod na- 
ture of the bandpass, we still determine a 
similar level of care to be necessary in the 
method. 
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